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AUDIO CODING SYSTEMS AND METHODS 

FIELD OF THE INVENTION 

This invention relates to audio coding systems and 
methods and in particular, but not exclusively , to such 
systems and methods for coding audio signals at low bit 
rates - 

BACKGROUND OF THE INVENTION 

In a wide range of applications it is desirable to 
provide a facility for the efficient storage of audio 
signals at a low bit rate so that they do not occupy large 
amounts of memory, for example in computers, portable 
dictation equipment, personal computer appliances, etc. 
Equally, where an audio signal is to be transmitted, for 
example to allow video conferencing, audio streaming, or 
telephone communication via the Internet, etc., a low bit 
rate is highly desirable. In both cases, however, high 
intelligibility and quality are important and this invention 
is concerned with a solution to the problem of providing 
coding at very low bit rates whilst preserving a high level 
of intelligibility and quality, and also of providing a 
coding system which operates well at low bit rates with both 
speech and music. 

In order to achieve a very low bit rate with speech 
signals it is generally recognised that a parametric coder 
or "vocoder" should be used rather than a waveform coder. 
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A vocoder encodes only parameters of the waveform, and not 
the waveform itself, and produces a signal that sounds like 
speech but with a potentially very different waveform. 

A typical example is the LPC 10 vocoder (Federal 
Standard 1015) as described in T.E. Tremaine "The Government 
Standard Linear Predictive Coding Algorithm: LPC10; Speech 
Technology, pp 40-4?, 1982) superseded by a similar 
algorithm LPClOe, the contents of both of which are 
incorporated herein by reference. LPC10 and other vocoders 
have historically operated in the telephony bandwidth (0- 
4kHz) as this bandwidth is thought to contain all the 
information necessary to make speech intelligible. However 
we have found that the quality and intelligibility of speech 
coded at bit rates as low as 2.4Kbit/s in this way is not 
adequate for many current commercial applications. 

The problem is that to improve the quality, more 
parameters are needed in the speech model , but encoding 
these extra parameters means fewer bits are available for 
the existing parameters. Various enhancements to the LPClOe 
model have been proposed for example in A.V. McCree and T.P. 
Barnwell III "A Mixed Excitation LPC Vocoder Model for Low 
Bit Rate Speech Coding"; IEEE-Trans Speech and Audio 
Processing Vol.3 No .4 July 1995, but even with all these the 
quality is barely adequate. 

In an attempt to further enhance the model we looked at 
encoding a wider bandwidth (0-8kHz) . This has never been 
considered for vocoders because the extra bits needed to 
encode the upper band would appear to vastly outweigh any 
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benefit in encoding it. Wideband encoding is normally only 
considered for good quality coders, where it is used to add 
greater naturalness to the speech rather than to increase 
intelligibility, and requires a lot of extra bits. 

One common way of implementing a wideband system is to 
split the signal into lower and upper sub-bands, to allow 
the upper sub-band to be encoded with fewer bits. The two 
bands are decoded separately and then added together as 
described in the ITU Standard G722 (X. Mai tre , K 7kHz audio 
coding within 64 kbit/s" , IEEE Journal on Selected Areas in 
Comm., vol.6, No. 2, pp283-298, Feb 1988). Applying this 
approach to a vocoder suggested that the upper band should 
be analysed with a lower order LPC than the lower band (we 
found second order adequate) . We found it needed a separate 
energy value, but no pitch and voicing decision, as the ones 
from the lower band can be used. Unfortunately the 
recombination of the two synthesized bands produced 
artifacts which we deduced were caused by phase mismatch 
between the two bands. We overcame 'this problem in the 
decoder by combining the LPC and energy parameters of each 
band to produce a single, high-order wideband filter, and 
driving this with a wideband excitation signal . 

Surprisingly, the intelligibility of the wideband LPC 
vocoder for clean speech was significantly higher compared 
to the telephone bandwidth version at the same bit rate, 
producing a DRT score (as described in W.D. Voiers, 
'Diagnostic evaluation of speech intelligibility' , in Speech 
Intelligibility and Speaker Recognition (M.E. Hawley, cd.) 
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pp. 374-387, Dowden, Hutchinson & Ross, Inc., 1977) of 86.8 
as opposed to 84.4 for the narrowband coder. 

However, for speech with even a small amount of 
background noise, the synthesised signal sounded buzzy and 
contained artifacts in the upper band. Our analysis showed 
that this was because the encoded upper band energy was 
being boosted by the background noise, which during the 
synthesis of voiced speech boosted the upper-band harmonics, 
creating a buzzy effect . 

On further detailed investigation we found that the 
increase in intelligibility was mainly a result of better 
encoding of the unvoiced fricatives and plosives, not the 
voiced sections . This led us to a different approach in the 
decoding of the upper band, where we synthesized only noise, 
restricting the harmonics of the voiced speech to the lower 
band only. This removed the buzz, but could instead add 
hiss if the encoded upper band energy was high, because of 
upper band harmonics in the input signal. This could be 
overcome by using the voicing decision, but we found the 
most reliable way was to divide the upper band input signal 
into noise and harmonic (periodic) components, and encode 
only the energy of the noise component. 

This approach has two unexpected benefits, which 
greatly enhance the power of the technique. Firstly, as the 
upper band contains only noise there are no longer problems 
matching the phase of the upper and lower bands, which means 
that they can be synthesized completely separately even for 
a vocoder. In fact the coder for the lower band can be 
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totally separate, and even be an off-the-shelf component. 
Secondly, the upper band encoding is no longer speech 
specific, as any signal can be broken down into noise and 
harmonic components, and can benefit from reproduction of 
the noise component where otherwise that frequency band 
would not be reproduced at all. This is particularly true 
for rock music, which has a strong percussive element to it. 

The system is a fundamentally different approach to 
other wideband extension techniques, which are based on 
waveform encoding as in McElroy et al : Wideband Speech 
Coding in 7.2KB/s ICASSP 93 pp 11-620 - 11-623. The problem 
of waveform encoding is that it either requires a large 
number of bits as in G722 (Supra) , or else poorly reproduces 
the upper band signal (McElroy et al) , adding a lot of 
quantisation noise to the harmonic components. 

In this specification, the term "vocoder" is used 
broadly to define a speech coder which codes selected model 
parameters and in which there is no explicit coding of the 
residual waveform, and the term includes coders such as 
mult i -band excitation coders (MBE) in which the coding is 
done by splitting the speech spectrum into a number of bands 
and extracting a basic set of parameters for each band. 

The term vocoder analysis is used to describe a process 
which determines vocoder coefficients including at least LPC 
coefficients and an energy value. In addition, for a lower 
sub-band the vocoder coefficients may also include a voicing 
decision and for voiced speech a pitch value. 
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SUMMARY OF THE INVENTION 

According to one aspect of this invention there is 
provided an audio coding system for encoding and decoding an 
audio signal , said system including an encoder and a 
decoder, said encoder comprising :- 

means for decomposing said audio signal into an upper 
and a lower sub* band signal; 

lower sub-band coding means for encoding said lower 
sub-band signal ; 

upper sub-band coding means for encoding at least the 
non-periodic component of said upper sub-band signal 
according to a source-filter model; 

said decoder means comprising means for decoding said 
encoded lower sub-band signal and said encoded upper sub- 
band signal , and for reconstructing therefrom an audio 
output signal , 

wherein said decoding means comprises filter means , and 
excitation means for generating an excitation signal for 
being passed by said filter means to produce a synthesised 
audio signal, said excitation means being operable to 
generate an excitation signal which includes a substantial 
component of synthesised noise in a frequency band 
corresponding to the upper sub-band of said audio signal . 

Although the decoder means may comprise a single 
decoding means covering both the upper and lower sub-bands 
of the encoder, it is preferred for the decoder means to 
comprise lower sub -band decoding means and upper sub- band 
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decoding means, for receiving and decoding the encoded lower 
and upper sub-band signals respectively. 

In a particular preferred embodiment, said upper 
frequency band of said excitation signal substantially 
wholly comprises a synthesised noise signal, although in 
other embodiments the excitation signal may comprise a 
mixture of a synthesised noise component and a further 
component corresponding to one or more harmonics of said 
lower sub-band audio signal. 

Conveniently, the upper sub-band coding means comprises 
means for analysing and encoding said upper sub-band signal 
to obtain an upper sub-band energy or gain value and one or 
more upper sub-band spectral parameters . The one or more 
upper sub-band spectral parameters preferably comprise 
second order LPC coefficients. 

Preferably, said encoder means includes means for 
measuring the noise energy in said upper sub-band thereby to 
deduce said upper sub-band energy or gain value. 
Alternatively, said encoder means may include means for 
measuring the whole energy in said upper sub-band signal 
thereby to deduce said upper sub-band energy or gain value. 

To save unnecessary usage of the bit rate, the system 
preferably includes means for monitoring said energy in said 
upper sub-band signal and for comparing this with a 
threshold derived from at least one of the upper and lower 
sub-band energies, and for causing said upper sub-band 
encoding means to provide a minimum code output if said 
monitored energy is below said threshold. 
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In arrangements intended primarily for speech coding , 
said lower sub-band coding means may comprise a speech 
coder, including means for providing a voicing decision. In 
these cases, said decoder means may include means responsive 
5 to the energy in said upper band encoded signal and said 
voicing decision to adjust the noise energy in said 
excitation signal dependent on whether the audio signal is 
voiced or unvoiced . 

Where the system is intended primarily for music, said 
10 lower sub-band coding means may comprise any of a number of 
suitable waveform coders, for example an MPEG audio coder. 

The division between the upper and lower sub-bands may 
be selected according to the particular requirements , thus 
it may be about 2.75kHz, about 4kHz, about 5.5kHz, etc. 
15 Said upper sub-band coding means preferably encodes 

said noise component with a very low bit rate of less than 
800 bps and preferably of about 300 bps. 

Where the upper sub-band is analysed to obtain an 
energy gain value and one or more spectral parameters, said 
20 upper sub-band signal is preferably analysed with relatively 
long frame periods to determine said spectral parameters and 
with relatively short frame periods to determine said energy 
or gain value . 

In another aspect, the invention provides a system and 
25 associated method for very low bit rate coding in which the 
input signal is split into sub-bands, respective vocoder 
coefficients obtained and then together recombined to an LPC 
filter. 
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Accordingly in this aspect, the invention provides a 
vocoder system for compressing a signal at a bit rate of 
less than 4.8Kbit/s and for resynthesizing said signal, said 
system comprising encoder means and decoder means, said 
encoder means including : - 

filter means for decomposing said speech signal into 
lower and upper sub-bands together defining a bandwidth of 
at least 5.5 kHz; 

lower sub-band vocoder analysis means for performing a 
relatively high order vocoder analysis on said lower sub- 
band to obtain vocoder coefficients representative of said 
lower sub-band; 

upper sub-band vocoder analysis means for performing a 
relatively low order vocoder analysis on said upper sub-band 
to obtain vocoder coefficients representative of said upper 
sub-band ; 

coding means for coding vocoder parameters including 
said lower and upper sub-band coefficients to provide a 
compressed signal for storage and/or transmission, and 

said decoder means including 

decoding means for decoding said compressed signal to 
obtain vocoder parameters including said lower and upper 
sub-band vocoder coefficients; 

synthesising means for constructing an LPC filter from 
the vocoder parameters for said upper and lower sub-bands 
and re- synthesising said speech signal from said filter and 
from an excitation signal. 

Preferably said lower sub-band analysis means applies 
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tenth order LPC analysis and said upper sub-band analysis 
means applies second order LPC analysis. 

The invention also extends to audio encoders and audio 
decoders for use with the above systems, and to 
corresponding methods . 

Whilst the invention has been described above it 
extends to any inventive combination of the features set out 
above or in the following description. 
BRIEF DESCRIPTION OF THE DRAWINGS 

The invention may be performed in various ways, and, by 
way of example only, two embodiments and various 
modifications thereof will now be described in detail" 
reference being made to the accompanying drawings, in 
which : - 

Figure 1 is a block diagram of an encoder of a first 
embodiment of a wideband codec in accordance with 
this invention; 

Figure 2 is a block diagram of a decoder of the first 
embodiment of a wideband codec in accordance with 
this invention; 

Figure 3 are spectra showing the result of the encoding- 
decoding process implemented in the first 
embodiment; 

Figure 4 is a spectrogram of a male speaker ; 

Figure 5 is a block diagram of the speech model assumed by 

a typical vocoder ; 
Figure 6 is a block diagram of an encoder of a second 

embodiment of a codec in accordance with this 
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invention; 

Figure 7 shows two sub-band short -time spectra for an 

unvoiced speech frame sampled at 16 kHz; 
Figure 8 shows two sub-band LPC spectra for the unvoiced 
5 speech frame of Figure 7; 

Figure 9 shows the combined LPC spectrum for the unvoiced 

speech frame of Figures 7 and 8; 
Figure 10 is a block diagram of a decoder of the second 
embodiment of a codec in accordance with this 
invention; 

Figure 11 is a block diagram of an LPC parameter coding 
scheme used in the second embodiment of this 
invention, and 

Figure 12 shows a preferred weighting scheme for the LSP 
predictor employed in the second embodiment of 
t hi s invent i on . 

In this description we describe two different 
embodiments of the invention, both of which utilise sub-band 
coding. In the first embodiment, a coding scheme is 
implemented in which only the noise component of the upper 
band is encoded and resynthesized in the decoder. 

The second embodiment employs an LPC vocoder scheme for 
both the lower and upper sub-bands to obtain parameters 
which are combined to produce a combined set of LPC 
parameters for controlling an all pole filter. 

By way of introduction to the first embodiment, current 
audio and speech coders, if given an input signal with an 
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extended bandwidth, simply bandlimit the input signal before 
coding. The technology described here allows the extended 
bandwidth to be encoded at a bit rate insignificant compared 
to the main coder. It does not attempt to fully reproduce 
the upper sub-band, but still provides an encoding that 
considerably enhances the quality (and intelligibility for 
speech) of the main bandlimited signal. 

The upper band is modelled in the usual way as an all- 
pole filter driven by an excitation signal. Only one or two 
parameters are needed to describe the spectrum. The 
excitation signal is considered to be a combination of white 
noise and periodic components, the latter possibly having 
very complex relationships to one another (true for most 
music) . In the most general form of the codec described 
below, the periodic components are effectively discarded. 
All that is transmitted is the estimated energy of the noise 
component and the spectral parameters; at the decoder, white 
noise alone is used to drive the all -pole filter. 

The key and original concept is that the encoding of 
the upper band is completely parametric - no attempt is 
made to encode the excitation signal itself. The only 
parameters encoded are the spectral parameters and an energy 
parameter. 

This aspect of the invention may be implemented either 
as a new form of coder or as a wideband extension to an 
existing coder. Such an existing coder may be supplied by a 
third party, or perhaps is already available on the same 
system (eg ACM codecs in Windows95/NT) . In this sense it 
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acts as a parasite to that codec, using it to do the 
encoding of the main signal, but producing a better quality- 
signal than the narrowband codec can by itself. An important 
characteristic of using only white noise to synthesize the 
5 upper band is that it is trivial to add together the two 
bands - they only have to be aligned to within a few 
milliseconds, and there are no phase continuity issues to 
solve. Indeed, we have produced numerous demonstrations 
using different codecs and had no difficulty aligning the 
10 signals. 

The invention may be used in two ways. One is to 
improve the quality of an existing narrowband (4kHz) coder 
by extending the input bandwidth, with a very small increase 
in bit rate. The other is to produce a lower bit rate coder 

15 by operating the lower band coder on a smaller input 
bandwidth (typically 2.75kHz) , and then extending it to make 
up for the lost bandwidth (typically to 5.5kHz) . 

Figures 1 and 2 illustrate an encoder 10 and decoder 12 
respectively for a first embodiment of the codec. Referring 

20 initially to Figure 1, the input audio signal passes to a 
low-pass filter 14 where it is low pass filtered to form a 
lower sub-band signal and decimated, and also to a high-pass 
filter 16 where it is high pass filtered to form an upper 
sub-band signal and decimated. 

25 The filters need to have both a sharp cutoff and good 

stop-band attenuation. To achieve this, either 73 tap FIR 
filters or 8th order elliptic filters are used, depending on 
which can run faster on the processor used. The stopband 
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attenuation should be at least 40dB and preferably 60dB, and 
the pass band ripple small - 0.2dB at most. The 3dB point 
for the filters should be the target split point (4kHz 
typically) . 

5 The lower sub-band signal is supplied to a narrowband 

encoder 18. The narrowband encoder may be a vocoder or a 
waveband encoder. The upper sub-band signal is supplied to 
an upper sub -band analyser 20 which analyses the spectrum of 
the upper sub-band to determine parametric coefficients and 
10 its noise component, as to be described below. 

The spectral parameters and the log of the noise energy 
value are quantised, subtracted from their previous values 
(i.e. differentially encoded) and supplied to a Rice coder 
22 for coding and then combined with the coded output from 
15 the narrowband encoder 18. 

In the decoder 12, the spectral parameters are obtained 
from the coded data and applied to a spectral shape filter 
23. The filter 23 is excited by a synthetic white noise 
signal to produce a synthesized non-harmonic upper sub-band 
20 signal whose gain is adjusted in accordance with the noise 
energy value at 24 . The synthesised signal then passes to 
a processor 26 which interpolates the signal and reflects it 
to the upper sub-band. The encoded data representing the 
lower sub-band signal passes to a narrowband decoder 30 
which decodes the lower sub-band signal which is 
interpolated at 32 and then recombined at 34 to form the 
synthesized output signal. 

In the above embodiment, Rice coding is only 
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appropriate if the storage/transmission mechanism can 
support variable bit-rate coding/ or tolerate a large enough 
latency to allow the data to be blocked into fixed- si zed 
packets. Otherwise a conventional quantisation scheme can be 
used without affecting the bit rate too much. 

The result of the whole encoding -decoding process is 
illustrated in the spectra in Figure 3, where the upper one 
is a frame containing both noise and strong harmonic 
components from Naklta by Elton John, and the lower one is 
the same frame with the 4 -8kHz region encoded using the 
wideband extension described above. 

Referring now in more detail to the spectral and noise 
component analysis of the upper sub-band, the spectral 
analysis derives two LPC coefficients using the standard 
autocorrelation method, which is guaranteed to produce a 
stable filter. For quantisation, the LPC coefficients are 
converted into reflection coefficients and quantised with 
nine levels each. These LPC coefficients are then used to 
inverse filter the waveform to produce a whitened signal for 
the noise component analysis. 

The noise component analysis can be done in a number of 
ways. For instance the upper sub-band may be full -wave 
rectified, smoothed and analysed for periodicity as 
described in McCree et al. However, the measurement is more 
easily made by direct measurement in the frequency domain. 
Accordingly, in the present embodiment a 256 -point FFT is 
perf ormed on the whitened upper sub-band signal . The noise 
component energy is taken to be the median of the FFT bin 
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energies. This parameter has the important property that if 
the signal is completely noise, the expected value of the 
median is just the energy of the signal. But if the signal 
has periodic components, then so long as the average spacing 
is greater than twice the frequency resolution of the FFT, 
the median will fall between the peaks in the spectrum. But 
if the spacing is very tight, the ear will notice little 
difference if white noise is used instead. 

For speech (and some audio signals) , it is necessary to 
perform the noise energy calculation over, a shorter interval 
than the LPC analysis. This is because of the sharp attack 
on plosives, and because unvoiced spectra do not move very 
quickly. In this case, the ratio of the median to the energy 
of the FFT, i.e. the fractional noise component, is 
measured. This is then used to scale all the measured energy 
values for that analysis period. 

The noise/periodic distinction is an imperfect one, and 
the noise component analysis itself is imperfect. To allow 
for this, the upper sub-band analyser 20 may scale the 
energy in the upper band by a fixed factor of about 50%. 
Comparing the original signal with the decoded extended 
signal sounds as if the treble control is turned down 
somewhat. But the difference is negligible compared to the 
complete removal of the treble in the unextended decoded 
signal . 

The noise component is not usually worth reproducing 
when it is small compared to the harmonic energy in the 
upper band, or very small compared to the energy in the 
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lower band. In the first case it is in any case hard to 
measure the noise component accurately because of the signal 
leakage between FFT bins. To some degree this is also true 
in the second case because of the finite attenuation in the 
5 stopband of the low-band filter. So in a modification of 
this embodiment the upper sub-band analyser 20 may compare 
the measured upper sub-band noise energy against a threshold 
derived from at least one of the upper and lower sub-band 
energies and, if it is below the threshold, the noise floor 

10 energy value is transmitted instead. The noise floor energy 
is an estimate of the background noise level in the upper 
band and would normally be set equal to the lowest upper 
band energy measured since the start of the output signal . 
Turning now to the performance of this embodiment, 

15 Figure 4, is a spectrogram of a male speaker. The vertical 
axis, frequency, stretches to 8000Hz, twice the range of 
standard telephony coders (4kHz) . The darkness of the plot 
indicates signal strength at that frequency. The horizontal 
axis is time. 

20 It will be seen that above 4kHz the signal is mostly 

noise from fricatives or plosives, or not there at all. In 
this case the wideband extension produces an almost perfect 
reproduction of the upper band. 

For some female and children's voices, the frequency at 

25 which the voiced speech has lost most of its energy is 
higher than 4kHz. Ideally in this case, the band split 
should be done a little higher (5.5kHz would be a good 
choice) . But even if this is not done, the quality is still 
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better than an unextended codec during unvoiced speech, and 
for voiced speech it is exactly the same. Also the gain in 
intelligibility comes through good reproduction of the 
fricatives and plosives, not through better reproduction of 
the vowels, so the split point affects only the quality, not 
the intelligibility. 

For reproduction of music, the effectiveness of the 
wideband extension depends somewhat on the kind of music. 
For rock/pop where the most noticeable upper band components 
are from the percussion, or from the ^softness" of the voice 
(particularly for females) , the noise-only synthesis works 
very well, even enhancing the sound in places. Other music 
has only harmonic components in the upper band - piano for 
instance. In this case nothing is reproduced in the upper 
band. However, subjectively the lack of higher frequencies 
seems less important for sounds where there are a lot of 
lower frequency harmonics. 

Referring now to the second embodiment of the codec 
which will be described with reference to Figures 5 to 12 
this embodiment is based on the same principles as the well- 
known LPCIO vocoder (as described in T. E. Tremain "The 
Government Standard Linear Predictive Coding Algorithm: 
LPCIO"; Speech Technology, pp 40-49, 1982), and the speech 
model assumed by the LPCIO vocoder is shown in Figure 5. 
The vocal tract, which is modeled as an all -pole filter 110, 
is driven by a periodic excitation signal 112 for voiced 
speech and random white noise 114 for unvoiced speech. 

The vocoder consists of two parts, the encoder 116 and 
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the decoder 118. The encoder 116, shown in Figure 6, splits 
the input speech into frames equally spaced in time- . Each 
frame is then split into bands corresponding to the 0-4 kHz 
and 4-8 kHz regions of the spectrum. This is achieved in a 
5 computationally efficient manner using 8th-order elliptic 
filters. High-pass and low-pass filters 12 0 and 122 
respectively are applied and the resulting signals decimated 
to form the two sub-bands. The upper sub-band contains a 
mirrored form of the 4-8 kHz spectrum. Ten Linear 

10 Prediction Coding (LPC) coefficients are computed at 124 
from the lower sub-band, and two LPC coefficients are 
computed at 126 from the high-band, as well as a gain value 
for each band. Figures 7 and 8 show the two sub-band 

short-term spectxa and the two sub-band LPC spectra 

15 respectively for a typical unvoiced signal at a sample rate 
of 16 kHz and Figure 9 shows the combined LPC spectrum. A 
voicing decision 128 and pitch value 130 for voiced frames 
are also computed from the lower sub-band. (The voicing 
decision can optionally use upper sub-band information as 

20 well) . The ten low-band LPC parameters are transformed to 
Line Spectral Pairs (LSPs) at 132, and then all the 
parameters are coded using a predictive quantiser 134 to 
give the low-bit -rate data stream. 

The decoder 118 shown in Figure 10 decodes the 

25 parameters at 13 6 and, during voiced speech, interpolates 
between parameters of adjacent frames at the start of each 
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pitch period. The ten lower sub-band LSPs are then 
converted to LPC coefficients at* 138 before combining them 
at 140 with the two upper sub-band coefficients to produce 
a set of eighteen LPC coefficients. This is done using an 
Autocorrelation Domain Combination technique or a Power 
Spectral Domain Combination technique to be described below. 
The LPC parameters control an all -pole filter 142, which is 
excited with either white noise or an impulse-like waveform 
periodic at the pitch period from an excitation signal 
generator 144 to emulate the model shown in Figure 5. 
Details of the voiced excitation signal are given below. 

The particular implementation of the second embodiment 
of the vocoder will now be described. For a more detailed 
discussion of various aspects, attention is directed to L. 
Rabiner and R.W. Schafer, 'Digital Processing of Speech 
Signals', Prentice Hall, 1978, the contents of which are 
incorporated herein by reference. 
IiPC Anal vs is 

A standard autocorrelation method "is used to derive the 
LPC coefficients and gain for both the lower and upper sub- 
bands. This is a simple approach which is guaranteed to 
give a stable all -pole filter; however, it has a tendency to 
over-estimate formant bandwidths. This problem is overcome 
in the decoder by adaptive formant enhancement as described 
in A.V. McCree and T.P. Barnwell III, 'A mixed excitation 
Ipc vocoder model for low bit rate speech encoding' , IEEE 
Trans. Speech and Audio Processing, vol.3, pp. 242-250, July 
1995, which enhances the spectrum around the formants by 
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filtering the excitation sequence with a bandwidth- expanded 
version of the LPC synthesis (all -pole) filter. To reduce 
the resulting spectral tilt, a weaker all -zero filter is 
also applied. The overall filter has a transfer function 
5 H(z) =A(z/0.5) /A(z/0 .8} , where A{z) is the transfer function 
of the all -pole filter. 
Resvnthesis LPC Model 

To avoid potential problems due to discontinuity 
between the power spectra of the two sub-band LPC models , 
10 and also due to the discontinuity of the phase response, a 
single high- order resynthesis LPC model is generated from 
the sub-band models . Prom this model , for which an order of 
18 was found to be suitable, speech can be synthesised as in 
a standard LPC vocoder. Two approaches are described here, 
the second being the computationally simpler method. 

In the following, subscripts L and H will be used to 
denote features of hypothesised low-pass filtered versions 
of the wide band signal respectively, (assuming filters 
having cut-offs at 4 kHz, with unity* response inside the 
pass band and zero outside), and subscripts 1 and h used to 
denote features of the lower and upper sub-band signals 
respectively. 

Power Spectral Domain Combination 

The power spectral densities of filtered wide-band 
signals P L (u) and P„(o)) , may be calculated as: 




' =1 a/(h)e-J wn | 2 . ifw<7r 

if 7T < < 2tt ; 



(1) 
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and 

where (n) , <\^(n) and gr,, gr,, are the LPC parameters and gain 
respectively from a frame of speech and p Jt p h , are the LPC 
model orders . The term tt-u/2 occurs because the upper sub- 
band spectrum is mirrored. 

The power spectral density of the wide-band signal, 
P w {o)) , is given by 

JVM = P L {u)-+Ph{v)- (3) 

The autocorrelation of the wide-band signal is given by 
the inverse discrete-time Fourier transf orm of P„(o>) , and 
from this the (18th order) LPC model corresponding to a 
frame of the wide-band signal can be calculated. For a 
practical implement at ion , the inverse transform is performed 
using an inverse discrete Fourier transform (DFT) . However 
this leads to the problem that a large number of spectral 
values are needed (typically 512) to give adequate frequency 
resolution , resulting in excessive computational 
requirements. 

Autocorrelation Domain Combination 

For this approach, instead of calculating the power 
spectral densities of "low-pass and high-pass versions of the 
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wide -band signal, the autocorrelations, r^{r) and r H (r) , are 
generated. The low-pass filtered wide-band signal is 
equivalent to the lower sub-band up- sampled by a factor of 
2. In the time-domain this up-sampling consists of 
5 inserting alternate zeros (interpolating) , followed by a 
low-pass filtering. Therefore in the autocorrelation 
domain, up- sampling involves interpolation followed by 
filtering by the autocorrelation of the low-pass filter 
impulse response. 

10 The autocorrelations of the two sub-band signals can be 

efficiently calculated from the sub-band LPC models (see for 
example R.A. Roberts and C.T. Mullis, 'Digital Signal 
Processing', chapter 11, p. 527, Adda son -Wesley , 1987). If 
rj (m) denotes the autocorrelation of the lower sub-band, then 

15 the interpolated autocorrelation, r'j(m) is given by: 

rl(m) = / r <(™/2> if m = 0, ±2, ±4, . . . (4) 
I 0 otherwise. 

The autocorrelation of the low-pass filtered signal r L (m) , 

is: 

r L (m) = r[{Tn) * (h(m) * h(—m)), (5) 



20 



where h(m) is the low-pass filter impulse response. The 
autocorrelation of the high-pass filtered signal r„(m) , is 
found similarly, except that a high-pass filter is applied. 

The autocorrelation of the wide-band signal r„(m) , can 
be expressed: 
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r w {m) = r L {m) + r H (m), (6) 

and hence the wide-band LPC model calculated. Figure 5 
shows the resulting LPC spectrum for the frame of unvoiced 
speech considered above. 
5 Compared with combination in the power spectral domain, 

this approach has the advantage of being computationally 
simpler. FIR filters of order 30 were found to be 
sufficient to perform the upsampling. In this case, the 
poor frequency resolution implied by the lower order filters 
10 is adequate because this simply results in spectral leakage 
at the crossover between the two sub-bands . The approaches 
both result in speech perceptually very similar to that 
obtained by using an high- order analysis model on the wide- 
band speech. 

15 From the plots for a frame of unvoiced speech shown in 

Figures 7 , 8, and 9, the effect of including the upper-band 
spectral information is particularly evident here, as most 
of the signal energy is contained within this region of the 
spectrum . 

20 Pitch/Voicing Analysis 

Pitch is determined using a standard pitch tracker. 
For each frame determined to be voiced, a pitch function, 
which is expected to have a minimum at the pitch period, is 
calculated over a range of time intervals. Three different 

25 functions have been implemented, based on autocorrelation, 
the Averaged Magnitude Difference Function (AMDF) and the 
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negative Cepstrum. They all perform well; the most 
-computationally efficient function to use depends on the 
architecture of the coder's processor. Over each sequence 
of one or more voiced frames, the minima of the pitch 
5 function are selected as the pitch candidates. The sequence 
of pitch candidates which minimizes a cost function is 
selected as the estimated pitch contour. The cost function 
is the weighted sum of the pitch function and changes in 
pitch along the path. The best path may be found in a 

10 computationally efficient manner using dynamic programming. 

The purpose of the voicing classifier is to determine 
whether each frame of speech has been generated as the 
result of an impulse-excited or noise-excited model. There 
is a wide range of methods which can be used to make a 

15 voicing decision. The method adopted in this embodiment 
uses a linear discriminant function applied to; the low-band 
energy, the first autocorrelation coefficient of the low 
(and optionally high) band and the cost value from the pitch 
analysis. For the voicing decision to work well in high 

20 levels of background noise, a noise tracker (as described 
for example in A. Varga and K. Ponting, 'Control Experiments 
on Noise Compensation in Hidden Markov Model based 
Continuous Word Recognition', pp. 167-170, Eurospeech 89) can 
be used to calculate the probability of noise, which is then 

25 included in the linear discriminant function. 
Parameter En coding 
Voicing Decision 

The voicing decision is simply encoded at one bit per 
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frame. It is possible to reduce this by taking into account 
the correlation between successive voicing decisions, but 
the reduction in bit rate is small. 
Pitch 

For unvoiced frames, no pitch information is coded. 
For voiced frames, the pitch is first transformed to the log 
domain and scaled by a constant (e.g. 20) to give a 
perceptually- acceptable resolution. The difference between 
transformed pitch at the current and previous voiced frames 
is rounded to the nearest integer and then encoded. 
Gains 

The method of coding the log pitch is also applied to 
the log gain, appropriate scaling factors being 1 and 0.7 
for the low and high band respectively. 
IiPC Coefficients 

The IjPC coefficients generate the majority of the 
encoded data. The LPC coefficients are first converted to 
a representation which can withstand quantisation, i.e. one 
with guaranteed stability and low * distortion of the 
underlying formant frequencies and bandwidths. The upper 
sub-band LPC coefficients are coded as reflection 
coefficients, and the lower sub-band LPC coefficients are 
converted to Line Spectral Pairs (LSPs) as described in F. 
JtaJcura, 'Line spectrum representation of linear predictor 
coefficients of speech signals', J". Acoust. 5oc. Ameri., 
vol.57, S35(A), 1975. The upper sub-band coefficients are 
coded in exactly the same way as the log pitch and log gain, 
i.e. encoding the difference between consecutive values, an 



WO 98/52187 



PCT/GB98/01414 



27 

appropriate scaling factor being 5.0. The coding of the 
low-band coefficients is described below. 
Rice Coding 

In this particular embodiment, parameters are quantised 
5 with a fixed step size and then encoded using lossless 
coding. The method of coding is a Rice code (as described 
in R.F. Rice & J.R. Plaunt, 'Adaptive variable -length coding 
for efficient compression of spacecraft television data', 
IEEE Transactions on Communication Technology, vol.19, 
10 no. 6, pp. 889-897, 1971), which assumes a Laplacian density of 
the differences. This code assigns a number of bits which 
increases with the magnitude of the difference. This method 
is suitable for applications which do not require a fixed 
number of bits to be generated per frame, but a fixed bit- 
15 rate scheme similar to the LPClOe scheme could be used. 
Voiced Excitation 

The voiced excitation is a mixed excitation signal 
consisting of noise and periodic components added together. 
The periodic component is the impulse* response of a pulse 
20 dispersion filter (as described in McCree et al) passed 
through a periodic weighting filter. The noise component is 
random noise passed through a noise weighting filter. 

The periodic weighting filter is a 2 0th order Finite 
Impulse Response (FIR) filter, designed with breakpoints (in 
2 5 kHz) and amplitudes: 
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The noise weighting filter is a 20th order FIR filter 
with the opposite response, so that together they produce a 
uniform response over the whole frequency band. 

LPC Parameter Encoding 

In this embodiment prediction is used for the encoding 
of the Line Spectral pair Frequencies (LSFs) and the 
prediction may be adaptive. Although vector quantisation 
could be used, scalar encoding has been used to save both 
computation and storage. Figure 11 shows the overall coding 
scheme. In the LPC parameter encoder 146 the input l s (t) is 
applied to an adder 148 together with the negative of an 
estimate l i (t) from the predictor 150 to provide a prediction 
error which is quantised by a quantiser 152 . The quantised 
prediction error is Rice encoded at 154 to provide an 
output, and is also supplied to an adder 156 together with 
the output from the predictor 150 to provide the input to 
the predictor 150. 

In the LPC parameter decoder 158, the error signal is 
Rice decoded at 160 and supplied to an adder 162 together 
with the output from a predictor 164. The sum from the 
adder 162, corresponding to an estimate of the current LSF 
component, is output and also supplied to the input of the 
predictor 164 . 

LSF Prediction 

The prediction stage estimates the current LSF 
component from data currently available to the decoder. The 
variance of the prediction error is expected to be lower 
than that of the original values, and hence it should be 
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possible to encode this at a lower bit rate for a given 
average error. 

Let the LSF element i at time t be denoted l^t) and the 
LSF element recovered by the decoder denoted l^it) . If the 
5 LSFs are encoded sequentially in time and in order of 
increasing index within a given time frame , then to predict 
l i (t) # the following values are available: 

OiWIl < j < i} 

and 

10 ih( T )\ T < ^ and 1 < j < 10}. 

Therefore a general linear LSF predictor can be written 
t~l 10 i_i 

where a i:J (T) is the weighting associated with the prediction 
of li(t) from X,(t-T) . 

15 In general only a small set of values of a^(T) should 

be used, as a high-order predictor is computationally less 
efficient both to apply and to estimate- Experiments were 
performed on unquantized LSF vectors (i.e. predicting from 
ljir) rather than X,(t) , to examine the performance of 

2 0 various predictor configurations, the results of which are: 
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Table 1 

System D (shown in Figure 12) was selected as giving the 
best compromise between efficiency and error. 

A scheme was implemented where the predictor was 
adaptively modified. The adaptive update is performed 
according to: 

C£ +I > =(l-/,)Ci*»+/«ixT 
C<* + "=(l-p)CW+p 3 /,x„ (8) 

where p determines the rate of adaption (a value of p=0.005 
was found suitable, giving a time constant of 4.5 seconds). 
The terms and are initialised from training data as 



and 




Here y ± is a value to be predicted (I^t)) and Xi is a vector 
of predictor inputs (containing 1, l^t-l) etc.). The 
updates defined in Equation (B) are applied after each 
frame, and periodically new Minimum Mean-Squared Error 
(MMSE) predictor coefficients ,p, are calculated by solving 
C J «P=C jey . 

The adaptive predictor is only needed if there are 
large differences between training and operating conditions 
caused for example by speaker variations, channel 
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differences or background noise. 
r>,iaTifcia ation aTt ^ Coding 

Given a predictor output l 4 (t) , the prediction error is 
calculated as e^ti-l^t) -li<t) . This is uniformly quantised 
by scaling to give an error i,(t) which is then losslessly 
encoded in the same way as all the other parameters. A 
suitable scaling factor is 160.0. Coarser quantisation can 
be used for frames classified as unvoiced. 
Results 

Diagnostic Rhyme Tests (DRTs) (as described in W.D. 
Voiers, 'Diagnostic evaluation of speech intelligibility', 
in Speech Intelligibility and Speaker Recognition (M.E. 
Hawley, cd.) pp. 374-387, Dowden, Hutchinson & Ross, Inc., 
1977) were performed to compare the intelligibility of a 
wide-band LPC vocoder using the autocorrelation domain 
combination method with that of a 4800 bps CELP coder 
(Federal Standard 1016) (operating on narrow-band speech) - 
For the LPC vocoder, the level of quantisation and frame 
period were set to give an average bit rate of approximately 
2400 bps. From the results shown in Table 2, it can be seen 
that the DRT score for the wideband LPC vocoder exceeds that 
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CELP 
Wideband LPC 


83.8 
86.8 



Table 2 

This second embodiment described above incorporates two 
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recent enhancements to LPC vocoders, namely a puise 
dispersion filter and adaptive spectral enhancement, but it 
is en^hasised that the embodiments of this invention may 
incorporate other features from the many enhancements 
published recently. 



